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The Ginzburg-Landau- Wilson theory that describes the 
disordered metal - d-wave superconductor phase transition 
at zero temperature is derived at weak coupling. The theory 
represents an interacting dissipative system of bosonic Cooper 
pairs in an effective random potential. I show that there ex- 
ists a wide crossover regime in the theory controlled by a line 
of Gaussian fixed points, each of which in two dimensions is 
characterized by a different universal value of the dc criti- 
cal conductivity. Relation to experiments on overdoped and 
underdoped cuprates is discussed. 

From the point of view of quantum phase transitions 
0, 1^, 1^, the BCS superconducting transition at zero 
temperature (T — 0) represents quite an unusual ex- 
ample. While the thermal superconducting transition 
in principle falls into the universality class of the XY- 
model, at r = Fermi liquid suffers the BCS insta- 
bility at an infinitesimal attractive interaction, which 
completely eliminates pair fluctuations. As long as the 
non-superconducting state remains a metal, the quan- 
tum phase transition into an s-wave superconductor is 
governed by a Gaussian fixed point, even in the pres- 
ence of weak non-magnetic impurities and Coulomb in- 
teraction 1^. This peculiar situation is caused by the 
same infrared singularity responsible for the Cooper phe- 
nomenon, which here renders the T — Q action a non- 
local function of the superconducting order parameter 
(OP), and its Landau expansion ill-behaved. For an s- 
wave superconductor with weak disorder Anderson the- 
orem guarantees the same singularity still to appear , 
and the nature of the Fermi liquid - superconductor 
transition essentially remains unchanged by the presence 
of randomness. 

The situation may be radically different if the pertur- 
bations are pair-breaking, like in case of a d-wave super- 
conductor with non-magnetic impurities. The effect of 
disorder on a d-wave superconductor becomes a partic- 
ularly pressing issue in the context of high-temperature 
superconductors, where at T = the superconducting 
transition can be tuned by a change of doping or of 
the level of disorder in the system Such quan- 

tum phase transitions (presumably of superconductor- 
insulator, or of superconductor-metal type) would cor- 
respond to the T = end points of the superconduct- 
ing region in the standard doping-temperature phase di- 
agram of a high- Tc material In this Letter I obtain 
the effective action that describes the Fermi liquid - d- 
wave superconductor T = critical point in the limit of 
weak coupling. Such a theory could be expected to apply 



to the quantum superconductor- metal transition induced 
by overdoping a high-Tc material. My starting point 
is a well-known notion that, when disorder is present, 
the attractive interaction in the d-wave channel has to 
be finite to produce pairing even at T = I show 
that this enables one to derive a local T — Ginzburg- 
Landau- Wilson (GLW) action for the fluctuating d-wave 
OP, which represents a system of interacting dissipative 
bosons (Cooper pairs) in an effective single-particle ran- 
dom potential. The singularities that in the s-wave case 
yield a non-local T = GLW action Q turn out to be de- 
moted to irrelevant terms (in the renormalization group 
sense) by the d-wave symmetry of the superconducting 
state. In a general case without the particle-hole sym- 
metry the pair propagator has both a dissipative and 
a propagative term linear in frequency. I flnd that the 
quantum critical behavior in the GLW theory is governed 
by a line of Gaussian fixed points (all with same critical 
exponents z = 2 and v = 1/2), except in an exponen- 
tially narrow critical region. In the Gaussian regime the 
small propagative term in the pair propagator represents 
a marginal operator, and leads to a non-universal singu- 
lar part of the T = dc conductivity at the transition in 
two dimensions. I conclude by commenting on the exper- 
iments on underdoped and overdoped cuprates in light of 
these results. 

Consider the BCS pairing Hamiltonian for spin- 1/2 
fermions in a standard Gaussian single-particle random 
potential V{r), with V{f)V{r') =WS{f-r'), with the 
fermions interacting only via weak attractive potential in 
the d-wave channel (Coulomb repulsion neglected). For 
simplicity I focus on two dimensions (2D) but neglect the 
weak- localization effects [0 (see later) , so that the nor- 
mal state is assumed to be a disordered metal. A general- 
ization to dimensions d > 2 is then straightforward. Us- 
ing the standard replica trick |jll| , the disorder-averaged 
free energy at T = is then F = lim„^o(l — Z'^)/n-, with 
Z^ = J D[-^^,^ exp{~S), S = So + S„u + Sd^s, and: 

q.i.(7 

^ /■ dwdvdv' Vj^k, I tt , ^ , ^ 
= - j -J2^3 + q.v + A 

X'^l-A-k, -ly)■^^,-ai-k', -i^')*.,^(fc' + 9, ^' + A, (2) 

^ W f duduj' . 

Sdis = -^ 2^ j ^^^y S{ki + fca - fcs - ki) 

ki . . .k4^,i,j .(7 ,a' 

X -^l^, {ki, u;')*j-,^' (fc2, ^')*U(fc3, c^)*,,.(fci, o;), (3) 
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FIG. 1. Quadratic (a), replica-diagonal quartic (b), and 
replica-off-diagonal quartic (c) kernels in the GLW theory. 

with being Grassman variables. Index i = 1, ...n nu- 
merates replicas, a spin, and integrations over the contin- 
uous Matsubara frequencies are over the whole real axis. 
I work in units in which ?i = fc^ = 1, and assume = 
V cos29j: j:,, with V > (attraction), and 6j:j:, the angle 
between the two momenta, ^(g) = ((7^/2rn) — e/, where 
is the Fermi energy, and I neglect the dependence of the 
pairing interaction on amplitudes of the momenta. The 
action S represents a minimal model that contains both 
the d-wave superconducting instability and the random 
potential. 

The standard procedure j|] is to decouple Sint in the 
Cooper channel with a complex Hubbard-Stratonovich 
field 



A(fc, q, Lo) = V^^^A{q, uj) cos 26lr, 



(4) 



and then integrate out fermions, to be left with the func- 
tional integral over A(q,uj) only. If performed exactly, 
this would amount to exact rewriting of the partition 
function. Since this is practically impossible, it is useful 
to expand the resulting action in powers of the OP field 

5gw[A] =^if2(l)A*(l)A,(l)+ J2 ^4,.,(1,2,3) 



6\ 



:A*(1)A,(2)A*(3)A,(1 + 3 - 2) + 0(A<'), 



(5) 



where I shortened the notation {qi,uji} — > 1, and 
dropped the regular A = part of the action. The 
quartic kernel in principle has replica-diagonal and off- 



diagonal parts: K, 



K. 



diag \ 



K^^^ ■ The quantum- 
mechanical average of the fluctuating field A is propor- 
tional to the d-wave superconducting OP. With disorder 
the integration over fermions can not be performed ex- 
actly even for the uj = q = Q mode. The kernels K2, 



K4, ... in the GLW action S'GLiy[A], however, are deter- 
mined exclusively by the connected correlation functions 
of the disordered metal, and may be represented diagram- 
matically as in Fig. 1. Full lines stand for the fermion 
propagator dressed by the random potential, and wavy 
lines for the fully renormalized vertex Wren- I assume 
that the fermion Green's function is diagonal in replica 
indices and given by the standard expression for a disor- 
dered metal 



G-\q,uj) =iLj- ^{q) + (^/2T)sign(a;). 



(6) 



At a weak attraction, the phase transition will occur at 
weak disorder, so I take the scattering time in the self- 
consistent Born approximation, r = {2TTWAf)~^ , with 
TV = rn/27r being the constant density of states per spin 
in d = 2. Accordingly, the vertex is approximated by the 
sum of particle-particle ladders (Fig. 2), which at small 
frequencies and momenta yields: 



WreniQ,n,iy) 



w 



\2iy + n\T+{PQ^/2)' 



(7) 



The choice of ladder diagrams for the vertex ensures that 
for an s-wave superconductor one reproduces a divergent 
K2 oc -VJ\fln{lq/2) when q -> (at w = 0) 0. This 
guarantees that in the s-wave case the system at T = 
becomes unstable towards pair formation at an infinites- 
imal attraction, in accord with the Anderson theorem 
Within the ladder approximation, for the uniform 
OP mode the present calculation of K2 becomes equiv- 
alent to the standard mean-field approach To the 
lowest order in q and tu, the quadratic (Fig. la), and the 
replica-diagonal part of the quartic kernel (at zero exter- 
nal momenta and frequencies) (Fig. lb) in the T = 
action for the OP field are then 



K2{q,io) = l-gH2efT)+g{^)' + 



AefT 



u) + 0{uj\qHn{ql)), 



(0) 



5g2 



12 A/" 



(8) 



(9) 



where g — VAf/2, and to regularize the standard ultravi- 
olet divergence in K2 (0) I employed a high-energy cutoff 

« £/ ^ 1/r. I = VfT is the elastic mean-free path. 
The parameter S — (ej — Q)/ef measures the particle- 
hole asymmetry in the model. 

Few remarks on the derivation of the above results are 
in order. A great simplification that arises due to the d- 
wave symmetry of the OP is that in ladder approximation 
for the vertex the contributions of second diagrams in 
Figs, la and lb at zero external momenta vanish due to 
integrations over the angles of momenta in the loops. The 
leading frequency dependence of K2 in this approxima- 
tion is thus entirely determined by the particle-particle 
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bubble at Fig. la. At g 9^ 0, however, the insertion 
of the particle-particle ladder into the second diagram 
at Fig. la gives a non-analytic momentum dependence 
oc q'^hi^ql) (at uj = 0). This is, however, irrelevant at 
long-distances compared to the already present and reg- 
ular term. This again is a consequence of the as- 
sumed d-wave symmetry of the OP. The reader should 
recall that in case of an s-wave superconductor exactly 
the same ladder approximation for the vertex leads to a 
diverging K2 when g — > 0, as noted earlier. Finally, in 
the calculation of i^4(0) the third diagram in Fig. lb is 
finite and actually of the same order in W as the first, 
and is retained in the Eq. 9. Within the ladder approxi- 
mation all the higher-order diagrams for K^^^'^^iO) vanish 
exactly. 

An important new element brought by disorder is the 
replica-off-diagonal contribution to the quartic kernel 
given by the sum of goggles-like diagrams at Fig. Ic. 
Besides mixing replicas, it conserves frequency for each 
replica separately, similarly as in Sdis in the fermionic 
action. In our ladder approximation for Wren the second 
and all higher-order diagrams in Fig. Ic again vanish for 
zero external momenta, and one obtains: 

J^off-Mag^Q-^ = -2g^WT^F{kfl){l + 0{S/efT)), (10) 

where F{x) « (Stt/S) ln(x)/a;, for a; ^ 1 (weak-disorder 
limit). Rescaling the OP field to bring the coefficient of 
the |w| term in Eq. (8) to unity as $ = {grY^^A, and 
retaining only the leading terms in the gradient expan- 
sion (Eqs. (8), (9) and (10)), the GLW action for the 
fluctuating d-wave OP at T = in real space takes a 
suggestive form: 

^["^l^E / d'rdr[^nr,r)i\dr\ + {S/{4efT))dr)^,ir,r) 

i 

+ ^ I V<f . (f, r ) p - I (f, r) p -f A I (f , r ) 1 4] 

j rfVdTdT'|$,(f,r)p|$^-(f,r')P, (11) 

where by the operator \dr \ one means \uj\ in Matsubara 
space. This GLW action describes a system of bosons, 
with a mass-like parameter = 2r/l^, a "chemical po- 
tential" lib — (ln(2e/T) — g^'^)/T, and the delta-function 
repulsion of strength A = 5/(12A/'), coupled to a bath 
of gapless fermions. The last, replica-coupling, term in 
the Eq. (11) is precisely as if it derived from averag- 
ing of the boson free-energy with a single-particle ran- 
dom potential term / V{,(f^|$(r, t)P in the action, with 

Vb{r)Vb{r') = Wb5{r- f), where Wb = 4VFF(fc/0. 

The obtained form of the T = GLW action appears 
natural and likely to be independent of the details of the 
underlying electronic Hamiltonian. Finite scattering rate 
in the fermion Green's function (6) cuts off the infrared 
singularities in case of d-wave and regularizes the theory 
very much like the finite temperature would do [|2j . Thus 




FIG. 2. Particle-particle ladder for the vertex Wren- The 
sum of incoming momenta is Q, frequency of the upper lines 
is + fl, and of the lower lines —f. 

a local GLW theory. Besides the standard frequency and 
momentum dependent terms, the pair propagator con- 
tains also a dissipative term of the Caldeira-Leggett form, 
signifying that the pairs can decay into gapless fermions. 
Notice that the ratio of the coefficients of the dissipative 
{\uj\) and the propagative (iuj) terms is proportional to 
dimensionless conductivity of the normal phase, so that 
for a good metal e/r ^ 1 and the propagative term is 
negligible, similarly to the situation at finite-T at weak 
coupling ||T^. The repulsion between the pairs, A, is pro- 
portional to the average spacing between the electron's 
energy levels near e/, since its origin is in the Pauli prin- 
ciple. Finally, the bosonic pairs themselves see an effec- 
tive random potential, which derives from the underlying 
disorder for electrons. 

The critical behavior of the interacting theory (11) 
with both dissipation and disorder is not known, and 
poses a rather challenging problem. In the weak-coupling 
(or weak-disorder) limit, however, the effective disorder 
Wb seen by the bosons matters only in a narrow crit- 
ical region, and can be neglected in the first approxi- 
mation. To see this, let me first drop the propagative 
term in the action when e/r ^ 1. The singular part of 
the disorder-averaged free-energy of the bosonic system 
described by the action (11) depends then on two di- 
mensionless couplings: A — Xrab, and Wb = Wbm},/\fib\, 
where jib is the renormalized chemical potential for the 
pairs. Using the above results for the bosonic parameters, 
it follows that Wb becomes comparable to A only when 
\g — gc\/gc ~ exp(— 1/gc), where (?c ^ 1 is the critical 
coupling. At weak enough coupling one is thus provided 
with a wide crossover region where the last term in GLW 
(11) can be safely neglected. 

The system described by the GLW theory (11) with 
Wb = in general has two phases at T = 0: the con- 
densed phase, i. e. a d-wave superconductor, at /ifc > 0, 
and the normal phase at fib < 0, with the quantum 
transition from one phase into another corresponding 
to pair-breaking. I assumed the normal phase to be 
a disordered metal. Right in 2D the weak-localization 
corrections will turn the normal state into an Ander- 
son insulator at distances of the order of localization 
length Lioc oi lexp{kfl) In the replica formal- 

ism this effect can be included by coupling the OP to 
massless "diffusons" ||ll|, i. e. rotations of the elec- 
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tron self-energy in replica space which leave the action 
at A = invariant. At weak coupling, however, these 
effects can become significant only in an even narrower 
region \g — gc\l Qc oc exp(— exp(l/(7c)), i- e. when the 
superconducting correlation length becomes of order of 
Lioc- In the Gaussian crossover regime under consider- 
ation the weak-localization effects are thus completely 
negligible even in 2D. The only remaining dimensionless 
coupling A = 57r/(6e/T) ^ I is irrelevant in all d > 2, 
since power counting yields the dynamical critical expo- 
nent 2 = 2, and the system is above (or in 2D right at) its 
upper critical dimension. The correlation length critical 
exponents in the crossover region is therefore a Gaussian 
1/ = 1/2, and the usual logarithmic corrections to scal- 
ing should be expected in 2D. Assuming that the T = 
superconductor-metal transition in overdoped cuprates is 
described by the GLW theory (11), this, for example, im- 
plies that Tc oc (xc — xY'^ oc (x^ — x) for x k, Xc, where Xc 
is the critical value of doping. Similarly, the superfluid 
density should be ps oc {xc - a;)''(''+^"2) ^ (^^^ _ xf/'^, 

in d = 3. This gives ps (x. x^'^^-^^^'^ ^ in d = 3, which in 
principle is testable experimentally. 

Since the normal state here is assumed to be a metal, 
at T = the dc conductivity consists of a finite regular 
term and a singular contribution: a = (Jr + as, where at 
weak disorder cr,. — (2e/T)e^/ft.. Since the critical point 
at Wb = is Gaussian, one can explicitly compute the 
singular part of the T = dc conductivity using the GLW 
action (11) Right at the transition in 2D I find: 

where /(O) = (tt^ - 4)/(167r), or 0.117(2e)V/i |l|. For 

X « 1, fix) = /(0)(l-((12-^2)/(^2_4))^2+O(^4))^ 

and for x » 1, f{x) = l/(4a;) + 0{l/x'^). Remarkably, 
if the system is not particle-hole symmetric, the univer- 
sality of the singular part of the conductivity jl^ at the 
transition in 2D is lost because the propagative term rep- 
resents a marginal operator at the Gaussian fixed point. 
The function f{x) interpolates between the universal val- 
ues for the particle-hole symmetric theory [6 — 0) and for 
the theory without dissipation (e/r — > 0), when a* — 0. 
Although (12) is not the directly experimentally relevant 
conductivity, where typically one assumes zero-frequency 
and low temperature |l^, it is the simplest one to 
calculate and may be indicative of the general situation. 
One could expect a non-universal a* whenever there are 
both linear frequency terms present in the OP action. If 
such a theory happens to be describing a superconductor- 
insulator transition, the total a — ag could still remain 
non-universal at the transition, much like what happens 
in s-wave superconductors in magnetic field [Q. 

The effective disorder Wb in (11) is a relevant perturba- 
tion at the Gaussian fixed point for all d < 4. The critical 
region dominated by Wb also grows with increase of disor- 
der. At a larger disorder, if the trend found here persists, 
the relative importance of the propagative and the dissi- 



pative terms in the OP action could become reversed. For 
a bad metal, if the coefficient of the dissipative term be- 
comes small compared to the propagative term the GLW 
theory (11) would approach the action for dirty-bosons 
JlTf , [Q. It is not inconceivable that at strong disorder 
the underdoped cuprates may be describable by the GLW 
theory in the form of (11), except with a small dissipative 
term, which would make the phase transition closer to the 
superfluid-Bose glass universality class. This would ac- 
count naturally for the universal critical conductivity ob- 
served at the superconductor-insulator transition in Zn- 
substituted underdoped cuprates 0. 
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